Scaling laws for the decay of multiqubit entanglement 
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We investigate the decay of entanglement of generahzed A'^-particle Greenberger-Home-Zeihnger (GHZ) 
states interacting with independent reservoirs. Scaling laws for the decay of entanglement and for its finite-time 
extinction (sudden death) are derived for different types of reservoirs. The latter is found to increase with the 
number of particles. However, entanglement becomes arbitrarily small, and therefore useless as a resource, 
much before it completely disappears, around a time which is inversely proportional to the number of particles. 
We also show that the decay of multi-particle GHZ states can generate bound entangled states. 
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Introduction. Entanglement has been identified as a key 
resource for many potential practical applications, such as 
quantum computation, quantum teleportation and quantum 
cryptography [l]. Being it a resource, it is of fundamental 
importance to study the entanglement properties of quantum 
states in realistic situations, where the system unavoidably 
loses its coherence due to interactions with the environment. 
In this context a peculiar dynamical feature of entangled states 
has been experimentally confirmed for the case of two qubits 
(two-level systems) |i2J: even when the constituent parts of 
an entangled state decay asymptotically in time, entanglement 
may disappear at a finite time JS] H |6l H E ^ | . The phe- 
nomenon of finite-time disentanglement, also known as en- 
tanglement sudden death (ESD) ||2l|71[8l|9l, illustrates the fact 
that the global behavior of an entangled system, under the ef- 
fect of local environments, may be markedly different from 
the individual and local behavior of its constituents. 

Since the speed-up gained when using quantum-mechanical 
systems, instead of classical ones, to process information is 
only considerable in the limit of large-scale information pro- 
cessing, it is fundamental to understand the scaling properties 
of disentanglement for multiparticle systems. Important steps 
in this direction were given in Refs. IS] IH |5l. In particular, 
it was shown in Ref. |3| that balanced Greenberger-Horne- 
Zeilinger (GHZ) states, j^-) = (|0)^^ + |1)®^)/V^, subject 
to the action of individual depolarization lH], undergo ESD, 
that the last bipartitions to loose entanglement are the most 
balanced ones, and that the time at which such entanglement 
disappears grows with the number N of particles in the sys- 
tem. Soon afterwards it was shown in Ref. [51 that the first 
bipartitions to loose entanglement are the least balanced ones 
(one particle vs. the others), the time at which this happens 
decreasing with N. A natural question arises from these con- 
siderations: is the ESD time a truly physically-relevant quan- 
tity to assess the robustness of multi-particle entanglement? 

In this paper we show that, for an important family of 
genuine-multipartite entangled states, the answer is no. For 
several kinds of decoherence, we derive analytical expressions 
for the time of disappearance of bipartite entanglement, which 
is found to increase with A^. However, we show that the time 



at which bipartite entanglement becomes arbitrarily small de- 
creases with the number of particles, independently of ESD. 
This implies that for multi-particle systems, the amount of en- 
tanglement can become too small for any practical application 
long before it vanishes. In addition, for some specific cases, 
we characterize not only the sudden-death time of bipartite en- 
tanglement but we can also attest full separability of the states 
in question. As a byproduct we show that in several cases the 
action of the environment can naturally lead to bound entan- 
gled states 1 10 1, in the sense that, for a period of time, it is not 
possible to extract pure-state entanglement from the system 
through local operations and classical communication, even 
though the state is still entangled. 

The exemplary states we take to analyze the robustness of 
multipartite entanglement are generalized GHZ states: 

|vl/o)^a|0)<^^ + /3|l)«~, (1) 

with a and /5 e C such that [ap+l/Jp = 1. Therefore, 
our results also constitute a generalization of those of Refs. 

Although ([T| represents just a restricted class of states, 
the study of its entanglement properties is important in its own 
right: these can be seen as simple models of the Schrodinger- 
cat state ifTTl . they are crucial for communication problems 
fl2l, and such states have been experimentally produced in 
atomic and photonic systems of up to iV = 6 lfT3l . 

Decoherence models. We consider three paradigmatic 
types of noisy channels: depolarization, dephasing, and a 
thermal bath at arbitrary temperature (generalized amplitude- 
damping channel). We consider N qubits of ground state |0) 
and excited state 1 1) without mutual interaction, each one indi- 
vidually coupled to its own noisy environment. The dynamics 
of the i-th qubit, 1 < i < TV, is governed by a master equation 
that gives rise to a completely positive trace-preserving map 
(or channel) £i describing the evolution as pi = EiPOi, where 
and Pi are, respectively, the initial and evolved reduced 
states of the i-th subsystem. 

The generalized amplitude-damping channel (GAD) is 
given, in the Born-Markov approximation, via its Kraus rep- 
resentation as |[Tl|9l 

£f ^^p, = EoP.eI + E^p.eI + E2P^El + E^p.eI- (2) 
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with i?o ^ V^^''^^^'^' + yr^|l)(l|), EE 

y^^|0)(l|, E2 EE y^(Vr^|0)(0| + |1)(1|) and 

£"3 = ^ ^-^^-^ p\l){0\ being its Kraus operators. Here n is 
the mean number of excitations in the bath, p = p{t) = 
1 — e~2'''*^^"+^^* is the probability of the qubit exchanging a 
quantum with the bath at time t, and 7 is the zero-temperature 
dissipation rate. Channel (j2]) is a generalization to finite tem- 
perature of the purely dissipative amplitude damping chan- 
nel (AD), which is obtainen from Q in the zero- temperature 
limit n = 0. On the other hand, the purely diffusive case is 
obtained from (j2| in the composite limit n ^ cx), 7 — > 0, and 
nj = r, where T is the diffusion constant. 

The depolarizing channel (D) describes the situation in 
which the i-th qubit remains untouched with probability l—p, 
or is depolarized - meaning that its state is taken to the maxi- 
mally mixed state (white noise) - with probability p. It can be 
expressed as 

ffA. = (l-p)p. + (p)l/2, (3) 

where 1 is the identity operator. 

Finally, the phase damping (or dephasing) channel (PD) 
represents the situation in which there is loss of quantum 
information with probability p, but without any energy ex- 
change. It is defined as 

Sf^p, = (l-p)p,+p(|0)(0|p.|0)(0| + |l)(l|p,|l)(l|). (4) 

The parameter p in channels Q, ^ and Q is a convenient 
parametrization of time: p — refers to the initial time and 
p = I refers to the asymptotic t ^ 00 limit. 

The density matrix corresponding to state ([T]i, 
po EE \^o){^o\ ^ |ap(|0)(Or^ + + 
Q/3*(|0)(l|)®^ + a*/3(|l)(0|)^^, evolves in time into a 
mixed state p given simply by the composition of all N 
individual maps: p = £182 ■■■ £nPq, where, in what follows, 
£i will either be given by Eqs. (|2]l, Q or 

Entanglement sudden death. In order to pick up the en- 
tanglement features of the studied states we will use the neg- 
ativity as a quantifier of entanglement |14|, defined as the 
absolute value of the sum of the negative eigenvalues of the 
partially transposed density matrix. In general, the negativ- 
ity fails to quantify entanglement of some entangled states 
(those ones with positive partial transposition) in dimensions 
higher than six 1 15 1. However, for the states considered here, 
their partial transposes have at most one negative eigenvalue, 
and the task of calculating the negativity reduces to a four- 
dimensional problem. So, in the considered cases, the nega- 
tivity brings all the relevant information about the separability 
in bipartitions of the states, i.e., null negativity means separa- 
bility in the corresponding partition. 

Application of channel (|2]) to every qubit multiplies the off- 
diagonal elements of po by the factor (1 — p)^^^, whereas 
application of channels ^ or by the factor (1 — p)^. The 
diagonal terms (|0)(0|)®^ and in turn give rise to 



new diagonal terms of the form (|0) (Ol)*^^-*^ (g) 
for 1 < fc < A^, and all permutations thereof, with coefficients 
Xk given below. In what follows we present the main results 
concerning the entanglement behavior of these states. 

Generalized amplitude- damping channel: Consider a bi- 
partition k : N ~ k of the quantum state. For channel (j2|, 
the coefficients A^^-° are given by X'^^^ = \a\'^x'^~''y'^ + 

\P\2^N-k^k^ With < X ^ 2^ + 1, y ^ J^, ^ 

^^I^^j^l-* and z = "^l"^"*"-* + 1 < 1. From them, the mini- 
mal eigenvalue of the states' partial transposition, A'^^^{p), 
is immediately obtained for the generahzed amplitude damp- 
ing channel IIT6I : 

Af^^(p) ^ 4 - ^Sl Ak . (5) 

Here 4 = l/2[X'^^^{p) + A^^f(p)] and A, = 
X'i^^{p)X%i^{p) - |a/3|2(l -p)^. Fromjs) one can see 
that|AP^(p)| < |AG^^(p)| < ... < |Af^^(p)|, for AT 

even, and |Af^^(p)| < |AG^^(p)| < ... <|A^(p)|,for 
A^ odd. 

The condition for disappearance of bipartite entanglement, 
A'^^^{p) = 0, is a polynomial equation of degree 2N. In 
the purely dissipative case n = 0, a simple analytical solution 
yields the corresponding critical probability for the amplitude- 
decay channel, pf^ (with f3 ^ 0): 

p^^(fc)=min{l,|a//3p/^}. (6) 

For \a\ < \(3\ probability (j6|l is always smaller than 1, mean- 
ing that bipartite entanglement disappears before the steady 
state is asymptotically reached. Thus, contition (j6]) is the di- 
rect generalization to the multiqubit case of the BSD condition 
of Refs. II2II2I for two qubits subject to amplitude damping. A 
remarkable feature about contition (j6]) is that it displays no de- 
pendence on the number of qubits k of the sub-partition. That 
is, the negativities corresponding to bipartitions composed of 
different numbers of qubits all vanish at the same time, even 
though they follow different evolutions. In the appendix we 
prove that at this point the state is not only separable accord- 
ing to all of its bipartitions but it is indeed fully separable, i.e., 
it can be written as a convex combination of product states. 

For arbitrary temperature, it is enough to consider the case 
k ~ N/2, as the entanglement corresponding to the most bal- 
anced bipartitions is the last one to disappear (we take A^ even 
from now on just for simplicity). For arbitrary temperature, 
the condition A'^^^{p) = reduces to a polynomial equation 
of degree A^, which for the purely diffusive case yields: 

p^'ff{N/2) = 1 + 2|a/3|2/^ - \/l + 4|a/3|4/A^ . (7) 

Depolarizing channel: For channel (Bj), the coefficients 
of p are given by Af = \a\^{l - + - 

IfiD^-KOnc obtains again Af (p) = 4 - VW^k, 
with 4 = l/2[Af + A^_,] and A, = Af A^_, - |a/3p(l - 
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■ Also here it is easy to show that the negativity asso- 
ciated to the most balanced bipartition is always higher than 
the others, while the one corresponding to the least balanced 
partition is the smallest one. The critical probability for the 
disappearance of entanglement in the N/2 : N/2 partition is 
given by: 

pf(7V/2) = l-(l+4|a/3|2/^)-i/^ (8) 

Phase damping channel: Finally, for the phase damping 
channel, whereas the off-diagonal terms of the density matrix 
evolve as mentioned before, all the diagonal ones remain the 
same, with = = A^^^. for 1 < < TV. In this case, 
A^^{p) = — |a/3|(l — p)^ . This expression is independent 
of fc, and therefore of the bipartition, and for any a, (3 ^ 
it vanishes only for p = 1, i. e., only in the asymptotic 
time limit, when the state is completely separable: general- 
ized GHZ states of the fonn (|7J, subject to individual dephas- 
ing, never experience BSD. 

The environment as a creator of bound entanglement. 
Some effort has been recently done in order to understand 
whether bound entangled (i.e. undistillable) states naturally 
arise from natural physical processes [ 171. In this context, it 
has been found that different many-body models present ther- 
mal bound entangled states 1 17]. Here we show, in a concep- 
tually different approach, that bound entanglement can also 
appear in dynamical processes, namely decoherence. 

For all channels here considered, the property |Ai(p)| < 
|A2(p)| < ... < |Aiv(p)| holds. Therefore, when |Ai(p)| = 
0, there may still be entanglement in the global state for some 
time afterwards, as detected by other partitions. When this 
happens, the state, even though entangled, is separable accord- 
ing to every 1 : — 1 partition, and then no entanglement can 
be distilled by (single-particle) local operations. 

An example of this is shown in Fig. [T] where the negativity 
for partitions 1 : A^ — 1 and N/2 : N/2 is plotted versus p, 
for AT = 4 and a = 1 /V2 = /3, for channel D. After the 1 : 3 
negativity vanishes, the 2 : 2 negativity remains positive until 
p — p^ (2) given by Eq. ([8]). Between these two values of 
p, the state is bound entangled since it is not separable but no 
entanglement can be extracted from it locally. Therefore, the 
environment itself is a natural generator of bound entangle- 
ment. Of course, this is not the case for channels AD and PD, 
since for the former the state is fully separable at p^^{k) (see 
Eq. (|6| and Appendix) while the latter never induces ESD. 

Does the time of ESD really matter for large N? Inspec- 
tion of critical probabilities (j6]), (j7]i and ([8]) shows that in all 
three cases p^ grows with A^. In fact, in the limit N ^ oo we 
have, for \af3\ ^ 0, p^^(fc) -> 1, p^'ff{N/2) -> 3 - w 
0.76 and p^{N/2) ^ 1 - -i^ w 0.55. This might be in- 
terpreted as the state's entanglement becoming more robust 
when the system's size increases. However, what really mat- 
ters is not that the initial entanglement does not disappear but 
that a significant fraction of it remains, either to be directly 
used, or to be distilled without an excessively large overhead 
in resources. The idea is clearly illustrated in Fig.|2] where 
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Figure 1 : Negativity as a function of p for a balanced, a = 1/ \/2 = 
13, four-qubit GHZ state and independent depolarizing channels. A 
similar behavior is observed with channel GAD with n 7^ 0, but the 
effect is not so marked (the smaller n, the weaker the effect). 



the negativity corresponding to the most balanced partitions 
is plotted versus p for different values of A^. Even though the 
ESD time increases with A^, the time at which entanglement 
becomes arbitrarily small decreases with it. The channel used 
in Fig.[2]is the depolarizing channel, nevertheless the behavior 
is absolutely general, as discussed in the following. 

For an arbitrarily small real e > 0, and all states for which 
\aj3\ ^ 0, the critical probability p^ at which h.j^/2{P() — 
eA7v/2(0), becomes inversely proportional to A^ in the limit of 
large A^. For channel (j2]), this is shown by letting k — N/2 
in which simplifies to A-%ffip) = -\af3\{l - p)^^^ + 

|Q,|5^Af/2yAf/2 ^ |^|2^Ar/2^Ar/2^ p^j. ^^^^ ^^^^ ^^(jj gj^-^j. 




P 



Figure 2: Negativity versus p for TV = 4, 40 and 400, for channel 
D and for the most balanced partitions. In this graphic a — 1/3 and 
13 = V8/3, but the same behavior is displayed for all other parame- 
ters and maps. The inset shows a magnification of the region in which 
vanishes. Even though |A|J)(p)| and |A^)o(p)l cross the lat- 
ter and vanish much later, they become orders of magnitude smaller 
than their initial value long before reaching the crossing point. 
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tation n, x^/^ and z^^^ are at most of the same order of 
magnitude as (1 — p)^/"^, whereas y^^^ and w^/'^ are much 
smaller than one. Therefore, for all states such that |a/3| ^ 
we can neglect the last two terms and approximate Q, at 

k = N/2, as A%ff{p) = -|a/3p(l - p)^/^. We set now 

A^fi'iPe) = eA%fi^{0) ^ e = (1 - p,)^/^ ^ log(e) = 

Y log(l -_Pe). Sincere < p^'^^{N/2) < 1, we can approx- 
imate the logarithm on the right-hand side of the last equal- 
ity by its Taylor expansion up to first order in p^ and write 
log(e) = -fpe, implying that pf^^ w -(2/iV) log(e). 
Similar reasonings applied to channels ^ and Q lead to 
p^'^^{t) Ri — (1/7V) log(e). These expressions assess the ro- 
bustness of the state's entanglement better than the ESD time. 
Much before ESD, negativity becomes arbitrarily small. The 
same behavior is observed for all studied channels, and all co- 
efficients a, P 0, despite the fact that for some cases, like 
for instance for channel Q, no ESD is observed. The pres- 
ence of log e in the above expression shows that our result is 
quite insensitive to the actual value of e <C 1 . 

Conclusions. We probed the robustness of the entangle- 
ment of generalized GHZ states of arbitrary number of par- 
ticles, N, subject to independent environments. The states 
possess in general longer ESD time, the bigger N, but the 
time at which such entanglement becomes arbitrarily small is 
inversely proportional to N. The latter time characterizes bet- 
ter the robustness of the state's entanglement than the time at 
which ESD itself occurs. In several cases the action of the 
environment can naturally lead to bound entangled states. An 
open question still remains on how other genuinely multipar- 
tite entangled states, such as graph states, behave. W states 
are expected to be more robust, since they have always only 
one excitation, regardless of TV IIT^ . For example, it is possi- 
ble to show that, for W states, channel AD induces no ESD; 
however, the negativity of the least balanced partitions decays 
with I/Vn |19|. This is another instance in which the ESD 
time is irrelevant to assess the robustness of multi-particle en- 
tanglement. Our results suggest that maintaining a significant 
amount of multiqubit entanglement in macroscopic systems 
might be an even harder task than believed so far. 

We thank F. Mintert and A. Salles for helpful comments and 
FAPERJ, CAPES, CNPQ, Brazilian Millenium Institute for 
Quantum Information, EU QAP project, Spanish MEC under 
FIS2004-05639, and Consolider-Ingenio QOIT projects for fi- 
nancial support. 

Appendix. Here we prove that the amplitude damping 
channel leads the state ([T]) to a fully separable state when all 
of its bipartite entanglements vanish. 

The evolved state can be written as p = |ap(|0)(0|)^^ + 
Ps, where ps is an unnormalized state. The goal is to show 
that Ps is fully separable. This is done by showing that ps 
is obtained, with a certain probability, from a fully separable 
state (7 through a local positive-operator-valued measurement 
(POVM) IT]. Because only local operations are applied, we 
conclude that ps, and thus p, must be fully separable. 

The (unnormalized) state a is defined as a = 2~^|/3p{l + 



(|0)(1|)®^ + f;(|l)(0|)«~]}, being 1 the 2^ x 2~ 
identity matrix. State a is GHZ-diagonal (see definition in 
Ref [201) and all of its negativities are null, then <t is 
fully separable ||20ll . Consider, for each qubit i, the local 
POVM with elements ^ = S{^/pf{k)\0) (0| + 

where S is such that A^^^^J'^ < 1, and 

A^*'^^^*' = 1 - A'i^^ A['1 Applying this POVM to every 
qubit of state a yields ps when the measurement outcome is 
TO = 1 (corresponding to Ai) for every qubit. □ 
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